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CN ' Abstract 

o: 

C^^ . We discuss the viability of a light particle (~ 30 eV neutrino) with strong self- 

^^ I interactions as a dark matter candidate. The interaction prevents the neutrinos from 

i-G ■ free-streaming during the radiation dominated regime so galaxy sized density perturba- 

II . tions can survive. Smaller scale perturbations are damped due to neutrino diffusion. We 

O I calculate the power spectrum in the imperfect fluid approximation, and show that it is 

'^ ' damped at the length scale one would estimate due to neutrino diffusion. The strength 

Cd ■ of the neutrino-neutrino coupling is only weakly constrained by observations, and could 

be chosen by fitting the power spectrum to the observed amplitude of matter density 

perturbations. The main shortcoming of our model is that interacting neutrinos can not 

;2j ■ provide the dark matter in dwarf galaxies. 



1 introduction 

Observations on many scales imply that most of the matter in the universe is not emitting 
electromagnetic radiation |jl| . The ratio of the energy density in this unseen matter component 
(usually refered as dark matter or DM) to the critical density today, is defined as flo and is 
under debate. Various dynamical determinations of Qo on the scales of galaxies and clusters 
give i7o ~ .2 ± .1 []T|. Observations of large velocity flows suggest that Qq could be larger than 
this, possibly fio — 1 with large uncertainties [0, 0. Theoretical prejudice favours Qo = 1, 
because this is more natural that fio < 1, and because it is predicted by inflation. However, in 
this case, one needs ^2 ~ .8 smooth on galaxy scales. 

Primordial nucleosynthesis implies that the energy density in baryons satisfies .003 < Qb < 
.06 [0, which might just be consistent with the observations of Qo- However, it is very difficult 
to build models involving only baryons that produce the large scale structure and the Cosmic 
Microwave Background (CMB) temperature fluctuations we see today P, |§], so it is usually 
assumed that the dark matter is made of something other than baryons. The list of potential 
candidates is long; the popular possibilities are massive neutrinos, axions, and supersymmetric 
particles. 

In the standard picture, the dark matter (DM) is in thermal equilibrium with the rest of 
the Universe at some early epoch, but effectively non-interacting when length scales of interest 
for structure formation come into the horizon. In this case, the only free parameter is the 
mass of the DM particle. If it had m ^ 1 GeV, it would be non relativistic at all times 
of interest for structure formation, and it would behave like a pressureless fluid. If it was 
light, it would be relativistic until cluster scales came within the horizon, and perturbations on 
smaller scales would be damped by free-streaming. As a result, in this second type of model 
(called hot dark matter or HDM models) clusters are produced first and must fragment into 
galaxies. This "top-down" galaxy formation scenario is difficult to reconcile with observations 
which indicate that galaxies and small scale structure formed first. A model where the 
dark particles behave non-relativistically is therefore prefered. The standard cold dark matter 
(CDM) model was found to have too much power on small scales when normalized to the 
amplitude of the measured COBE/DMR temperature anisotropics and several modifications 
were introduced. Another possibility is a "warm" dark matter particle with a mass (~ keV) 
intermediate between that of "hot" (~ eV) and "cold" (> GeV) 0. Currently, models with a 
mixture of hot and cold dark matter (MDM) are favoured @] . 

It is usually assumed that during structure formation (T^ < 100 eV) the evolution of the 
Universe was controlled by electromagnetism and gravitation, while new particle physics takes 
place at accelerator energies. However, it is possible that the dark matter particles could have 
some unexpected interaction during the epoch of structure formation, despite this occuring at 
a low energy scale. Some examples of this have been previously studied [0, 11, 12, 1^. In this 



paper, we consider a simple interacting dark matter model; we assume that the dark matter 
consists of light particles with a large cross section for scattering off each other (we discuss later 
the meaning of "large"). For definiteness, we take the light particle to be a neutrino with mass 
~ 20 — 30 eV (as for HDM). Our interest is to determine whether, in the presence of hot dark 
matter self-interactions, density perturbations on scales smaller than galaxy clusters remain 
undamped. In normal HDM, they are washed out by free streaming of relativistic particles; in 
an interacting model, one would rather expect the perturbations to be damped on a smaller 
scale by diffusion (Silk damping). We will show that this is in fact the case, and compute on 
what scale they are damped. 



The question of whether "sticky neutrinos behave hke cold dark matter" has been previously 
addressed by Raffelt and Silk [jlO[. They argue that light (m ~ 20 eV) interacting neutrinos will 



dissipate small scale perturbations (< galaxy sized) by "Silk damping", but there will be no 
free streaming. This result was briefly questioned by Dicus et al. in their paper about neutrino 



interactions in supernovae []I4|, where they suggest that cosmological density perturbations in 
a neutrino fluid would dissipate at the speed of sound [|. It is well known that perturbations in 
a relativistic perfect fluid in a radiation dominated Universe neither grow nor are damped [^ . 
However, if the mean free path of the neutrinos is non-zero, they do not constitute a perfect 
fluid. There are contributions to the fluid viscosity and heat conduction coefficients that are 
proportional to the mean free path. The perturbations in the relativistic neutrino fluid in the 



early Universe will be damped by the fluid imperfections, and as Weinberg |T^ showed for Silk 
damping, the damping length scale can be estimated by the random walk argument used in [|l^ . 
Raffelt and Silk's result is therefore a correction to a perfect fluid analysis, not a contradiction 
to it. 

The aim of this paper is to study "sticky neutrinos" as dark matter in the imperfect fluid 
approximation. Providing that the particle mean free path is much shorter than the lengthscale 
of the density perturbations, this is a sensible approximation because it takes into account the 
non-zero mean free path of the particles, but avoids the complexity of the Boltzmann equation. 
In section 2, we briefly review the experimental bounds on the vv scattering cross section, and 
discuss theories in which such interactions could be strong. In the flrst part of the third section, 
we review the Boltzmann equation and its relation to the fluid approximation, in the interest 
of making this paper more self-contained. In the second part of this section, we estimate the 
scale on which perturbations in an "sticky hot dark matter" Universe would be damped. We do 
this by treating the neutrinos as a perfect fluid in an expanding Friedmann-Robertson- Walker 
Universe to lowest order, and then computing the viscosity as a perturbation proportional 
to the neutrino mean free path. We briefly discuss the subsequent behaviour of the sticky 
neutrinos in galaxy formation, and calculate the power spectrum. In section 4 we show that 
interacting neutrinos cannot constitute the halo of dwarf galaxies. However, we argue that this 
is a potential problem not only for light fermions but for all dissipationless DM. Finally, we 
present our conclusions. 

2 particle physics 

Despite the precision to which most Standard Model quantities are known, there are no ex- 
perimental bounds on the vv or vv scattering cross section. This is hardly surprising, as it 
would be difficult to scatter neutrino beams off each other with a high enough luminosity to 
see anything within the Standard Model. The best constraint on the vv and vv cross section 
comes from the observation of the neutrinos from SN1987A. These had to pass through the 
"cosmic neutrino background radiation" between here and the Large Magellanic Cloud (LMC) 
before arriving at the earth. Approximately the right number of neutrinos were detected, which 
places an upper bound on the vv interaction. This constraint has been computed in [O]; we 
briefly outline their arguments here. 

The "cosmic background" of light stable neutrinos has a number density rij, = ^{Ti,/T^)^n^, 
where T^/T^ depends on how many different families of particles annihilate into photons after 

""^They note, however, that their flat space results are not directly applicable to cosmology. 



the neutrino gas decouples from the photons. In this paper, we will assume {T^/T^Y = 4/11, 
n„ ~ n^/3. 

The flux of Ey ~ 10 MeV neutrinos from the supernova will arrive at the earth if the mean 
free path A of a supernova neutrino is of order the distance D ~ 2 x 10^^ cm from here to the 



LMC The mean free path is computed in [17] with some care, but can be roughly estimated as 

A-i ~ an, (1) 

where a is the cross section for a neutrino to scatter oS a u or a u. Requiring X/D ^ 1 gives 

a < 90 GeV"^ at v^ = 20 keV (2) 

where y'Ssn is the centre-of-mass energy of the collision. Since we assume that the neutrinos 
are the dark matter, they have masses of order 20 — 30 eV, and Ssn — ^Eyirti, ~ (20 keV)^. 

We are interested in the neutrino self interaction cross section at centre-of-mass energy 
scales near 10 eV. This is because galaxy scales come into the horizon when T ^ 10 eV, and 
the neutrinos become non-relativistic shortly thereafter; if galaxy scale perturbations are going 
to be damped, it will be during this period. We therefore need to scale @) down to lower 
values of s. iFiovn a phenomenological point of view, a can scale as a{s) = o'{ssn)s/ssn (due to 
exchange of a boson with m :^ s), or as a{s) = cr{ssn)ssn/s (exchange of a boson with m <^ s). 
Galaxy-sized perturbations will be undamped in both cases. We will concentrate on a model 
where the cross section scales as s. The other possibility, that it would scale as 1/s, would lead 
to an interaction strong enough to affect the dynamics of groups of galaxies today. Long range 
dark matter interactions have been considered in ||l^ . 



Sticky neutrino models are not easy to construct. The obstacle is that the Standard Model 
neutrinos are members of the same SU{2) doublet as the electrons, so it is difficult to give them 
strong self-interactions without contradicting present experimental data. This can be avoided 



by introducing right-handed gauge singlet neutrinos [0, but then one has difficulties with the 
primordial nucleosynthesis bound on the number of light degrees of freedom present at T ~ 1 
MeV in the early Universe 0]: Ni, ^ 4. Attempting to construct a neutrino mass matrix that is 
consistent with this constraint, and various indications of neutrino oscillations (solar neutrino 
problem, atmospheric neutrino deficit,...) requires more fine-tuning. 

Sticky neutrinos were originally studied in the context of a modified triplet majoron model 
P^ , PU[ , where the majoron (the nambu-goldstone boson associated with the spontaneous break- 
ing of lepton number L) was more massive than the lightest neutrino. The neutrinos interact 
via the exchange of the majoron. However, the triplet majoron couples to the Z, and has 
not been detected at LEP, so this model is ruled out. If one extends the Standard Model 
fermion content to include right-handed neutrinos, then these can acquire majorana masses 
by coupling to a singlet scalar with a vacuum expectation value, and Dirac masses with the 



left-handed neutrinos by coupling to the Higgs. This is the singlet majoron model [0], and is 
experimentally only weakly constrained. The angular component of the scalar is the majoron, 
and its coupling to the light neutrinos is proportional to the neutrino mass. This model can be 
fine-tuned to produce interacting light neutrinos. The right-handed neutrino majorana masses 
must be of order the Dirac masses, to allow a big enough neutrino scattering cross section by 
majoron exchange. To avoid having too many neutrinos at nucleosynthesis, one can make one 
family light enough that the right handed component is not yet in equilibrium, and another 
heavy enough to decay beforehand. (We have added strong neutrino self-interactions, so it is 
possible to arrange fast decays for a heavy neutrino). This adds up to three neutrinos plus the 



majoron (= .6z/) at nucleosynthesis. One must then impose flavour symmetries, to prevent the 
middle neutrino, who is the dark matter candidate, from decaying. This mass matrix can not 
explain the solar neutrino problem by oscillations. 

Alternatively, we can disconnect the exchanged scalar from the neutrino mass generation 
mechanism. In this way, we can get "three neutrinos at nucleosynthesis" . If we introduce a few 
MeV scalar with gauge strength couplings to the dark matter neutrino, it induces strong enough 
interactions to preserve galaxy scale perturbations, and may not be present as a relativistic 
degree of freedom at nucleosynthesis. We must again arrange the neutrino mass matrix so 
that the heaviest neutrino decays before nucleosynthesis while the middle one, who is our dark 
matter candidate, is stable, and the lightest one's right-handed component is not in thermal 
equilibrium at T ~ MeV. The light neutrino must therefore have a very small Dirac mass, and 
very weak coupling to the scalar. There is again no possibility of explaining the solar neutrino 
problem. 

We can construct a neutrino mass matrix that could explain the solar neutrino problem if we 
allow four neutrinos at nucleosynthesis. We introduce a ~ few MeV scalar whose interactions 
are unrelated to the neutrino masses, give it gauge strength couplings to the 30 eV neutrino, 
and very weak interactions with the other two neutrinos, who are lighter. We assume that 
these other two neutrinos also have very small Dirac masses, so their right-handed components 
are not in thermal equilibrium at nuclesynthesis. This makes four neutrinos at nucleosynthesis. 
We can then allow the left-handed components of the two lighter neutrinos to have a majorana 
mass difference that explains the solar neutrino deficit. 

Another possibility would be to disconnect the interacting hot dark matter from the neu- 
trinos, and imagine that there is a light real self-interacting scalar as well as three "standard" 
neutrinos. This would give 3.6 neutrino families at nucleosynthesis, and the neutrino mass 
matrix is still free, but one loses the great advantage of neutrinos over other dark matter candi- 
dates, which is that they are known to exist. In any case, we are mainly interested in interacting 
hot dark matter from a phenomenological point of view, and do not advocate any particular 
model. 

3 Evolution of density perturbations. 

The aim of this section is to discuss the formalism for calculating the evolution of linear per- 
turbations in interacting dark matter, and apply it to the case of "sticky neutrinos". In the 
first part, we review the Boltzmann equation and its relationship to the fluid approximation. 
We discuss how to compute the viscosity ?7, and the heat conduction x for a slightly imperfect 
fluid, and determine their contribution to the damping of perturbations. In the second part, 
we calculate rj for a neutrino fluid, and show that the estimate of the damping scale made by 



Raffelt and Silk |I0[ is essentially correct. 
We use the metric 

ds^ = dt^ — o?{t)[5ij + hij{x,t)]dx^dx^ . (3) 

where hij is the metric perturbation in synchronous gauge. We are interested in the effect 
of particle interactions on the evolution of matter density perturbations within the horizon. 
On those scales, metric perturbations are negligible compared to pressure gradients, damping 
and other hydrodynamical effects [^, |2^. We therefore neglect the metric perturbations in 
our analytic estimates, although they are included in our numerical calculation of the power 
spectrum. 



We define k to be the magnitude of tlie comoving tliree-niomentum, and 7^ its directional 
cosine: if \piP^\^^'^ = p is tlie magnitude of the physical momentum, p = k/a, pi = —k'ji, and 
p^ = a~^pji. Roman indices run from 1 to 3, and greek from to 3. 

3.1 The Boltzmann equation and the fluid approximation. 

The Boltzmann equation describes the evolution of the phase space density of the particles 
f{x\pi,t), and provides a framework in which one can calculate the growth and survival of 
density perturbations once they come within the horizon. It can be written as 

m = c[f] (4) 

where £[/] is the derivative of / along the path in phase space of the particles, and C[f] is 
the collision integral. If the particles described by / have no collisions, / is constant along 
a particle's path in phase space, or C[f] =0. In a relativistic fluid under the influence of 



gravitation, C[f] is defined as ||24 



^[/]=^^|^-re,i>V^. (5) 

C[f] parametrizes the changes in the distribution function due to the particle collisions. The 
integral representing the change in the distribution function fa for particle a due to the process 

a + b^c + d is ^ 



C[fa\ = \j (m,(mcdIiAMab^cd\\fah ' fcfd) (6) 

where |A1p is the matrix element squared for the process a + h ^ c + d, including spin sums 
and averages, and 



,5' = (2ir)*«*(p„+pi-p,- 

Implicitly, in equation (^ we have used Maxwell-Boltzmann statistics by neglecting the occu- 
pation in the final phase space. This approximation is justified if the chemical potentials of all 
particles are negligible at the temperature scales involved. But even in this limit, the solution 
of the Boltzmann equation is complicated. 

One can avoid dealing with the collision term by integrating (^ over quantities that are 
conserved in the interaction decribed by C[f] |^. This eliminates the right-hand-side of (^. 
For instance, if p^ is the four-momentum of the particle a, one can write 

J dUaP^^Llfa] = J dUaP^^Cifa] (9) 

and it is easy to show from symmetry considerations that the integrand on the right hand side 
is proportional to Pa + Pb ~ Pc ~ Pdi ^o the integral is zero. Therefore, one obtains 

dIiaP':,C[fa] = . (10) 



In the case of sticky neutrinos, the scattering interaction conserves particle number, so we also 
have 

ldUX[fa\=0. (11) 

Defining the stress-energy tensor for a single species as 



T^" = 2 dUfp^p" , (12) 

and the number density vector 

Nf" = 2 f dUfp^" , (13) 



it is easy to check, in a Friedmann-Robertson- Walker Universe, that ([TOD and (11) are equivalent 
to the usual conservation equations 

j.^u^ = Tf""^^ + T'^^T"^ + T^J'^'^ = (14) 

and 

N^,^ = Q. (15) 

It is clear, however, that since these are only conservation equations, they do not contain 
enough information to specify the dynamics. We still need to solve the Boltzmann equation 
d^) to obtain /. Since the Universe was homogeneous and isotropic to better than one part in 
10^ when it was radiation dominated, one can solve the Boltzmann equation for perturbations 



away from the homogenous and isotropic background ||T5| . A second possibility, which can only 
be used for interacting particles, is to perturb away from a perfect fiuid. We shall use this 
second approach for sticky neutrinos, since we assume they are strongly interacting. These 
two schemes are not identical, because in the first case the background is a homogeneous and 
isotropic equilibrium distribution / ~ exp{— £'/T(t)}, and the perturbation in / corresponds 
to a density perturbation. In the second case, one perturbs away from a locally thermal distri- 
bution / ~ exp{— i5/T(a;, t)}, so density perturbations are already present in the background 
solution. It is a two step procedure to find the evolution of the density perturbations in the fiuid 
approximation. We first estimate the magnitude of the perturbation away from local thermal 
equilibrium using the Boltzmann equation. Then we calculate the contributions of these per- 
turbations to the stress-energy tensor, and solve for the evolution of the density perturbations 
from the conservation equation ([l^)- The perturbations away from local thermal equilibrium 
introduce new terms in (0) which can affect the growth of the density perturbations. 

A perfect fluid is locally in thermal equilibrium, so we can write the phase space density for 
particles consituting an imperfect fluid as 



f = g^'\E,x\t)+g^^\E,p,,x\t) (16) 

where 

(27r)3 



g^''\E,x\t) = -^exp{-[E - fi{x,t)]/T{x,t)} (17) 



is the local equilibrium distribution, g is the number of degrees of freedom of the particle, and 
E = U"pa, where U" is the four-velocity of a particular observer, g'^^^ measures the departure 
of the distribution from exact local equilibrium. 



If the particles involved in an interaction are in local thermal equilibrium, then the collision 
term in the Boltzmann equation is zero; however, g^^^ does not neccessarily satisfy the left hand 
side of the Boltzmann equation. Assuming that g^^'^ <^ g^^\ we have 



j^dg^ ^ PJ^dg^ ^ ^^2d9^^^ 



dt 



dx' 



BE 



(niu(mjYidS^\M 



ab^cd 



A9^^+9f^^-9^^-9T^) 



(18) 
(19) 



In principle, to correctly compute the deviations from perfect fluid behaviour, we still have 
to evaluate the collision integral. However, if we only wish to estimate viscosity and heat 
conduction coefficients, the right hand side can be taken |^ (to within a factor of ~ 4) as 
its first term ~ EaVg'^^'^ where F is the thermally averaged interaction rate for the reaction 
described by C[f]. This gives 



9' 



ES 



(20) 



To get the damping scale of the perturbations correctly, we would need a better approximation, 
as can be found, for instance, in [Q. A more accurate approach would require a numerical 
integration of the coupled Einstein and Boltzmann equations that describe the behaviour of 
metric perturbations and the density field of the different particles [^ . 

If the mean time between particle collisions is r = T~^, and t is the timescale associated 
with the density perturbation in the fluid, ( pOf ) implies g^^^ ~ g^^^r /t. So g'^^^ -C g^^\ and the 
fluid approximation is consistent, provided that the size of the density perturbations in the 
early Universe is much longer than the mean free path of the particles. In this limit, the phase 
space density to flrst order in r is 



f = 9^'^ ^ ^C[g^'A 



(21) 



We have not yet deflned /i or T, or discussed in which reference frame we are measuring the 
energy E. In principle E = p°'Ua, so we must specify Ua- We follow Landau and Lifschitz [28], 



and take Ua to be the velocity associated with energy flux. (Note that this is not the choice 



made in [|16|, so we can not immediately compare to Weinberg's equations). This is natural in 
the approximation scheme we have used, because 



nOJ 



2 / dUEp'g'^^' = 2 / dUp'C[g(^>] = 0{t 



(0)1 







(22) 



where the last equality follows from (PUJ). We deflne T such that the equilibrium expression 
for the energy density is equal to the energy density in the reference frame where there is no 
energy flux, and /i such that the equilibrium expression for the number density is equal to the 
measured number density in the same reference frame. (In other words, the perturbations g^^^ 
will not contribute to p or n in the reference frame where T^^ = 0). 

The general form of the stress-energy tensor and the particle flux vector A^'^ for a slightly 
imperfect relativistic fluid can be determined by requiring the entropy to be non- decreasing 
In the rest frame of U, the velocity associated with energy flux, they are 
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-r/(9if/2 + d^U^) 


_^(5lf/3 + ^Sfjl) 
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-r]{d^U^ + d^U^) 


-7]{d^U^ + aif/3) 
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(23) 



and 



N^' 



n, 




d' 



(24) 



Formal expressions for the heat conduction coefficient x ^i^nd the bulk and ordinary viscosities ( 
and f] have been elegantly calculated in |^ for a generic relativistic fluid. However, their inte- 
grals are difficult to evaluate. We can simply approximate the fluid parameters by substituting 
(|2l| ) into ([T^) and ([13|) , and comparing to the formulae for T^'^ and A^'^. To extract x from A^^ 
one must use the constraint (|22D. We will do this for sticky neutrinos in the next subsection. 

Once T^'^ for an imperfect fluid has been calculated, the conservation equation ([I^) can be 
used to study the evolution of density perturbations on scales much longer than the particle 
interaction length. As noted in [|l^], perturbations in a fluid of relativistic particles are princi- 
pally damped by viscosity. If we neglect the metric perturbations (reasonable for sub-horizon 
sized modes) and the heat conduction coefficient, and decompose the density perturbation in 
Fourier modes, (|T^) implies 

S = -t9 (25) 



e + ^-e 

a 



-1^ 

p {- - 



P 



(26) 



where 6 = ikjU\ 6 = Sp/p, k is the comoving three-momentum and p = k/a is the physical 
three momentum. These are the equations for the longitudinal component { U \\ p) oi the 
perturbation; the transverse component is rapidly damped |jl6|. One can see from (^) and 
([26|) that the perturbation oscillates and damps at a rate 



2p'^rj 



Hp + p) 



(27) 



due to the fluid viscosity r]. Including heat conduction x and bulk viscosity ( gives ||16 

->2 



p- 



Hp + p) 

where A is a function of various fluid parameters 



C + -V + AX 



(2J 



A 



df 



-1 



■^— (i) -"=-(f 




(29) 



Vs denotes the speed of sound. This function A goes to zero as particles become extremely 
relativistic. 



3.2 Perturbation Damping for Sticky Neutrinos 

The length scale on which neutrino diffusion damps density perturbations increases with the age 
of the Universe until the neutrinos become non-relativistic. We can estimate the comoving scale 
on which density perturbations will be damped by treating the neutrinos as relativistic until 
T ~ mi,/3, and neglecting any further motion after this temperature. A similiar approximation 
provides the ordinary neutrino free streaming scale to within a factor of three, and since we only 
know the fluid viscosity and heat conduction to within the same order, it should be adequate. 



In the extreme relativistic limit, the bulk viscosity ( is zero, so we will neglect it. To 



compute the viscosity rj, we note from ||2^ that in the reference frame where U = (1, 0, 0, 0) 



\ dx2 dxi J 
Substituting (|21|) into (|T2]), we find 

r'^^^/p>.ny-,vf^i^ + ... (31) 

where the "..." indicate the terms which are zero when the angular integral is performed. With 

the identity 

/Att 
dn-fy-f^-f^ = — {5'^''^ + 5''' 5^^ + 5'^ 5^^) (32) 

15 

this gives in the relativistic limit 

4rnT , , 

V = ^ (33) 

which agrees with the result in |]16|, ^, ^ . 

The heat conduction coefficient x can be estimated by substituting (^) into (|13|), (with the 
constraint (|2^) ) and comparing to (plj) . In the relativistic limit, we get 



4™ , , , 



which is a factor of 3 smaller than the result in |T^, ^ but is within the errors of our approxi- 
mation. 

We can roughly estimate the physical length scale £d on which perturbations will be damped 
as in ^ 2tt/p£), where poit) is is determined from (^) by setting 

t-nr 

T{t)dt = 1 (35) 

and using the relativistic expression for rj (0). We neglect heat conduction since A ~ in 
the relativistic limit, and perturbations are principally damped by viscosity, tnr denotes the 
moment when neutrinos become non-relativistic and occurs when T ~ m^/?) ~ 10 eV. The 
physical lengthscale on which perturbations are damped is 

(36) 

This is approximately the distance one would estimate that a relativistic neutrino could travel 
via a random walk in time t: a neutrino can take N = t/r steps of length r, so on average 
each neutrino travels a distance yNr. One expects the neutrino number density to become 
homogeneous on scales shorter than ~ t/\/N, which was the estimated damping scale of neu- 
trino density perturbations given in |]TU[. To summarize, if one includes interactions for a hot 
dark matter candidate, perturbations are damped by diffusion rather than by free-streaming. 
The length scale on which perturbations are washed out is therefore shorter, and the power 




spectrum should look more like that of cold dark matter, as claimed in |10 . 

We have neglected mass effects in (|36|) , because our fluid parameters are only correct to 
within factors of 3 or 4. To properly determine the damping scale, we would need accurate 

10 



determinations of Xi V ci-nd (, valid as the temperature drops past the neutrino mass. We 
could then integrate F over all times to determine io- Instead, to check that our relativistic 
approximation is acceptable, we have calculated x ^iid t] in the non-relativistic limit, matched 
them to the relativistic versions at T = m^, and calculated ir,. The damping scale is (to within 
a factor of 2) the same as in (pBI). 

The comoving scale up to which density perturbations are washed out by neutrino diffusion 
is approximately 

XD = ^ = ^iD^ 3 /^ Mpc (37) 

kn Tq \j Mpc 

where r is the co-moving mean free path at T„r = 10 eV. As discussed in section 2, the neutrino 
mean free path is bounded below by the arrival of the neutrinos from SN1987A. This means 
that the scale on which perturbations are damped is also bounded below, and we need to check 
that galaxy-sized perturbations do indeed survive. 

The supernova bound implies that the timescale for neutrino self-interactions at tnr (when 
Tj, = 10 eV) must satisfy 



no CFsn i -LQ 



J n \ / Oxn \ 



Tnr ^ Tsn = [ T=^ ] [ T^n (38) 



where a subscript means the present value of a given magnitude, n is the neutrino number 
density, a is the cross section, and r^^ is the distance to the LMC. Snr and Sgn are respectively 
the centre-of-mass energy for the neutrino self-interaction at T^ = 10 eV, and for scattering 
a supernova neutrino off the neutrino CMB. The ratio Sgn/snr is to the power +1 if the cross 
section scales as s (exchanged boson mass > ^/s^), and to the power —1 for a cross section 
scaling as 1/s (exchanged boson mass < ^/s^). For a non-renormalisable cross section, the 
supernova bound implies that perturbations on comoving scales ^ 10~^ Mpc are damped by 
neutrino diffusion (or ~ 10"^ Mpc for a cross section that scales as 1/s). The supernova bound 
on the cross section is therefore consistent with the survival of galaxy scale perturbations, 
whether the neutrino stickiness is due to the exchange of light or heavy particles. 

3.3 Evolution in the non— relativistic regime and final power spec- 
trum. 

We have seen that the sticky neutrinos behave quite differently from standard neutrinos while 
they are relativistic. However, once they are non-relativistic, they become a pressureless dust 
to lowest order, like most other dark matter candidates. If the neutrino cross section scales as 
s, then the interactions become weaker as the temperature drops, and our fluid approximation 
breaks down. However, if the neutrino interaction timescale is of order the horizon size then 
the interaction is irrelevant. On the other hand, if our fluid approximation is still valid, then we 
know that the viscosity and heat conduction coefficients are proportional to the pressure [^, 
and the pressure is negligible. So we expect the non-relativistic, interacting neutrinos to behave 
like standard neutrinos, if the interaction is short range. A longer range interaction (exchanged 
boson mass -C m^) is potentially more complicated, particularily if the interactions are still 



strong today. Dark matter with long-range interactions has been considered in |]13[ ; we neglect 
it here. Constraints on dark matter self-interactions from galaxy formation and evolution in the 
non-linear regime have been estimated in [0. The authors argue that the average interaction 



length of dark matter today must be longer than 100 Mpc. Sticky neutrinos who interact 
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via the exchange of a m ~ MeV boson and who are consistent with the supernova bound are 
essentially non-interacting today (the interaction rate scales as a~^ or a~^), so they easily satisfy 
the constraints listed in [|12|. 

Based on our analytic results, we expect the power spectrum of interacting HDM to be sim- 
ilar to standard CDM on large and intermediate scales, but damped on small scales by neutrino 
diffusion. We are working in the fluid approximation, so to compute the power spectrum, we 
must numerically integrate the Einstein equations and the energy conservation equation ([I^ ) 



for an imperfect fluid. This is considerably simpler than integrating the coupled Einstein and 
Boltzmann equations. 

Let us briefly outline here the equations and approximations used. We perform our calcu- 
lations in the synchronous gauge using the metric (|]). The equations governing the evolution 
of perturbations are more easily written in Fourier space. If we only consider only scalar-type 
perturbations, the metric perturbations can be written in terms of two scalar fields h{k, t) and 



r]{k,t) as follows p7 



hij{x, t)= I £ke^^''[kikjh{k, t) + {hkj - -5ifiT]{k, t)}] . (39) 

As in the FRW case, the equations governing the evolution of matter density perturbations can 
be derived from the conservation equation (^1|). For a relativistic fluid, with viscosity but no 
bulk viscosity or heat conduction, equations ( pS]) and (p6D become: 

-5 = -Ue- ]-h] (40) 



3 V 2 

o + -e = p'{--\e-^)\ (41) 

a \^ P J 



This results differs from P| or P^ because our matter does not behave as a perfect fluid. As 



the temperature drops past the neutrino mass, the neutrinos no longer behave as relativistic 
particles. The evolution of density perturbations in the E r^ m regime is complicated by the 
time-dependent relation between the neutrino momentum and energy [^. This adds terms to 
(|40| ) and (^Tj) and makes the numerical calculation of the pressure and density more difficult 
p2| . To simplify the calculation, we assume the transition from a relativistic gas to presureless 
dust to be instantaneous. We also assume that neutrino interactions are frozen in the non- 
relativistic regime due to the small mean velocity. As a result, we treat the neutrino gas as 
CDM particles with zero velocity in the synchronous gauge, whose evolution is given by: 

5 = h (42) 

As previously noted, perturbations in a relativistic fluid are damped by viscosity. The off- 
diagonal spatial elements of the stress-energy tensor (parametrized by the viscosity) are also 
gauge-invariant, which makes them simple to calculate. We therefore include only the viscosity 
(no heat conduction coefficient and no bulk viscosity) in our numerical calculation of the power 
spectrum. We take the physical interaction time r to scale as a^ (comoving t ^ a^), and assume 
an initial Harrison-Zel'dovich spectrum of perturbations from inflation. 

In Fig.l we give the numerical integration of the previous set of differential equations coupled 
with the equations for the evolution of metric perturbations, baryons and photons. We used a 
modified version of the COSMIC package made publicly available by Berstchinger and Ma, and 
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described in [^ . The dashed hne corresponds to standard CDM and is plotted for comparison. 
The three sohd hnes correspond, in decreasing amphtude, to zero viscosity (zero mean free 
path) and comoving mean free paths at T^ = 10 eV of 10^^ Mpc and 0.1 Mpc. In the case of 
zero viscosity, the neutrino perturbations oscillate as acoustic waves around the CDM power 
spectrum. The behaviour is similar to the baryon-photon plasma where the interactions cause 
the perturbations in the baryon component to oscillate as sound waves. When the viscosity 
is introduced the oscillations are damped approximately at the scale given by the diffusion 
length. The effect of viscosity is to damp the power spectrum at the scale given by (^7\): 3 and 
7 Mpc— 1 respectively for co-moving mean free paths of 0.01 and 0.1 Mpc. The damping is a 
purely fluid effect and it has nothing to do with gravitation as was already remarked in |l^ . 

One should notice that the overall power spectrum amplitude is not completely accurate on 
scales k ~ 0.3Mpc~^. We assumed, when calculating the power spectrum, that the neutrino gas 
underwent an instantaneous transition from radiation to pressureless dust at T^ = 10 eV. For 
the purpose of showing the effect of damping on small scales due to viscosity, the instantaneous 
transition approximation is quite accurate, and simplifies the calculations. 

To summarize, neutrino viscosity helps to reduce power on small scales while keeping most 
of the features of CDM. In this way, it alleviates the problems that plague the standard CDM 
model, such as large scale pairwise velocity dispersion, that come from having too little power 
on large (> 30 h~^ Mpc) scales relative to small (< 10 h^^ Mpc) ones. A detailed analysis 
of the success of our model will require N-body simulations to compare with observations as 
performed, for example, by 0. However, such analysis goes beyond the scope of the present 
paper. 

4 Phase Space Constraints 

The main argument against neutrinos as DM candidates was put forward by Gunn and Tremaine 
PB| . They argued that the large phase space density of galaxies was close to the density of a 
degenerate Fermi gas and they virtually ruled out neutrinos as dark matter candidates. The 
argument has been refined since [^ ^ . Currently, the tightest constraints come from two 
dwarf spheroidal galaxies within our local group: Draco and Ursa Minor. In those galaxies, if 
the spatial distribution of dark matter and luminous stars is similar, the central dark matter 
densities needed to explain the kinematic data are poM = 0.5 — IMqPc^. For non-interacting 
neutrinos, the conservation of the phase space density implies [^ : 



,„„>irOe^fl^y'V^H2^y",;./. (43) 

\rDM J \ Pdm J 

To overcome this limit and bring rrii, down to 30 eV, the core radius of the DM, tdm, should 
be almost two orders of magnitude larger than the one of the luminous stars, making the mass 
to light ratio of the system close to 10^Mq/Lq. For sticky neutrinos the bound is very similar 
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However, the halo of dwarf galaxies is not only a difficult test for fermionic DM, but any 
dissipationless galaxy formation scenario. The density of an object that formed without dissi- 
pation uniquely specifies its redshift of formation 



1 + .,„ =. 4h-J" L^lJ (44) 
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Figure 1: Power spectrum for three different "sticky neutrino" models. The sohd hues corre- 
spond (in decreasing amphtude) to interactions with no neutrino diffusion, and with co-moving 
mean free paths at Ty = 10 eV of .01 Mpc and .1 Mpc. These give rise to estimated damping 
scales of 7 Mpc~^ and 3 Mpc~^, respectively. The dashed line correspond to standard CDM 
and is plotted for comparison. The y-axis scale is arbitrary. We took Ho = 50Km s— IMpc— 1. 



where pr, is the central density of the system. For Draco and Ursa Minor, the previous equation 
shows that (1 + Ztum) — 30. In CDM, MDM and similar scenarios [Q the first structures turn 
around at about (1 + ztum) — 10 and to bring down the previous limit would be necessary 
to have a DM core radius 4 times larger than the luminous matter and M/L ~ IOOOM0/L0 
for these dwarfs, which is also rather high. Lacking a clear understanding of star and galaxy 
formation one can not disregard the possibility that a significant fraction of the halo is baryonic 
or dissipational [^, which would weaken the limit on the neutrino mass. 



5 Discussion 

As mention in the introduction, the "standard" CDM model, when normalized to the amplitude 
of temperature anisotropics measured by COBE/DMR, has too much power on small scales. 
Several modifications were proposed that agree better with observations, such as, for instance, 
hot plus cold dark matter models, CDM models with cosmological constant, tilted power spec- 
trum and variants thereof. However, if one is going to add free parameters to one's dark matter 
model, it is in principle as reasonable to add interactions for a single component as to mix two 
(or three) components. As an example of this, we have considered a dark matter model consist- 
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ing of ~ 30 eV neutrino with a large vv scattering cross section. We have argued in this paper 
that these interacting neutrinos, although not easy to produce in a reasonable particle physics 
model, may be a viable dark matter candidate. Unlike ordinary hot dark matter, they do not 
free-stream. In section 3 we estimated the damping scale associated with neutrino diffusion 
and calculated the power spectrum. The current upper limits on the vv cross section permit 
a damping scale shorter than IMpc and accordingly galaxy scale density perturbations will 
survive. Stronger upper limits on the cross section will increase the damping scale and could 
invalidate the model. The main positive feature of our model is that the power on small scales 
is reduced by neutrino diffusion while on larger scales we retain the main feautures of standard 
CDM. The neutrino-neutrino cross section, which is the only free parameter, could possibly 
be tuned to give "the right amount of power on small scales". The main handicap is that 
interacting neutrinos, like ordinary neutrinos, can not constitute the halos of dwarf galaxies. 
However, it can be argued that dwarfs are difficult to explain in any dark matter scenario, so 
this may not be a fatal problem. 
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